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P2X receptors are trimeric cation channels that open in response to the binding of adenosine triphosphate (ATP) 
to a large extracellular domain. The x-ray structure of the P2X4 receptor from zebrafish (zfP2X4) receptor reveals 
that the extracellular vestibule above the gate opens to the outside through lateral fenestrations, providing a po- 
tential pathway for ions to enter and exit the pore. The extracellular region also contains a void at the central axis, 
providing a second potential pathway. To investigate the energetics of each potential ion permeation pathway, we 
calctilated the electrostatic free energy by solving the Poisson-Boltzmann equation along each of these pathways in 
the zfP2X4 crystal structure and a homology model of rat P2X2 (rP2X2). We found that the lateral fenestrations 
are energetically favorable for monovalent cations even in the closed-state structure, whereas the central pathway 
presents strong electrostatic barriers that would reqtiire structural rearrangements to allow for ion accessibility. 
To probe ion accessibility along these pathways in the rP2X2 receptor, we investigated the modification of introduced 
Cys residues by methanethiosulfonate (MTS) reagents and constrained structural changes by introducing disulfide 
bridges. Our results show that MTS reagents can permeate the lateral fenestrations, and that these become larger 
after ATP binding. Although relatively small MTS reagents can access residues in one of the vestibules within the 
central pathway, no reactive positions were identified in the upper region of this pathway, and disulfide bridges 
that constrain movements in that region do not prevent ion conduction. Collectively, these restilts suggest that ions 
access the pore using the lateral fenestrations, and that these breathe as the channel opens. The accessibility of 
ions to one of the chambers in the central pathway likely serves a regulatory function. 



INTRODUCTION 

P2X receptors are cation-permeable channels that play 
crucial roles in extracellular ATP-mediated signaling 
(North, 2002; Burnstock, 2006). These receptors are 
widely expressed throughout the human body and play 
important roles in cell-cell communication, ranging from 
sensory signal transduction (Cook et al., 1997; Finger 
et al., 2005) to immune responses (Solle et al., 2001; 
Chessell et al., 2005). Malfunction of P2X receptor- 
mediated ATP signaling is implicated in various patho- 
logical conditions, including neuropathic pain (Souslova 
et al., 2000; Dell'Antonio et al., 2002; Tsuda et al., 2003) 
and inflammatory diseases such as asthma and arthritis 
(Dowd et al., 1998; Idzko et al., 2007). Therefore, P2X 
receptor channels are promising targets for the devel- 
opment of new therapeutic agents. 

P2X receptor channels are composed of three sub- 
units that assemble as either homomeric or heteromeric 
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complexes of seven different subunits (P2X1-7) (Nicke 
et al., 1998; Vial et al., 2004) . Each of the three subunits 
of P2X receptors has two transmembrane helices, which 
link the intracellular amino and carboxyl termini to the 
large extracellular region (^^280 amino acids). These ar- 
chitectural features distinguish P2X receptors from the 
other ligand-gated cation channels, which include tet- 
rameric glutamate receptors (Rosenmund et al., 1998; 
Sobolevsky et al., 2009) and pentameric acetylcholine 
receptors (Reynolds and Karlin, 1978; Unwin, 2005). The 
recent crystal structure of a closed P2X4 receptor from 
zebrafish (zfP2X4) revealed that this large extracellular 
region has a novel architecture among the available ion 
channel structures, harboring three potential ATP-binding 
pockets, disulfide-rich motifs, and conspicuous vestibules 
that form avoid in the middle of three subunits (Kawate 
et al., 2009). Although accessibility and metal bridging 
experiments in rat P2X2 (rP2X2) support a TM2 helix- 
straightening model for opening of the pore (Kracun 
et al., 2010; Li et al., 2010), the accompanying structural 
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rearrangements in the extracellular region remain rela- 
tively unexplored. 

In this study, we investigated the pathway through 
which ions traverse the extracellular domain of P2X re- 
ceptors to enter or exit the transmembrane pore. The 
x-ray structure of zfP2X4 reveals that two potential path- 
ways exist. The first pathway is through three lateral fen- 
estrations that connect the extracellular solution to an 
extracellular vestibule just above the occluded gate region 
of the pore (Fig. 1, A and B). Although these lateral fen- 
estrations have dimensions of '^8 x 10 A in the zfP2X4 
receptor structure, several important side chains were 
not modeled because of the poor quality of electron 
density maps in this region (N54, D59, T60, Q329, 1335, 
and 1336 in zfP2X4) . The inclusion of all side chains in 
the zfP2X4 structure suggests that these fenestrations 
may be somewhat narrower. The second possible path- 
way for ions to access the pore is along the central three- 
fold axis of the zfP2X4 structure (Fig. 1 , A and B) . Two 
conspicuous electronegative cavities are present between 
the outermost end of the extracellular region and the 
extracellular vestibule in zfP2X4, which are predicted 
to form favorable sites for cation binding. Although the 
diameter of the pathways connecting these chambers with 
each other and the extracellular solution at the top of 
the receptor would not support rapid ion throughput, 
they might expand after the binding of ATP. There is 
precedence in other ion channels for both types of po- 
tential permeation pathways. Ions must traverse an ex- 
tended (^85 A) pore that projects into the cytoplasm in 
inward-rectifier potassium (Kir) channels (Nishida and 
MacKinnon, 2002; Nishida et al., 2007; Tao et al., 2009) , 
and ions are thought to use four lateral fenestrations 
between the transmembrane pore and the cytoplasmic 
"hanging gondola" in voltage-activated potassium (Kv) 
channels (Kobertz and Miller, 1999; Long et al., 2007), 
and in corresponding regions in CNG channels (Johnson 
and Zagotta, 2005). To probe the ion access pathway 
in the extracellular region of a P2X receptor, we per- 
formed Cys-scanning mutagenesis and accessibility studies 
using methanethiosulfonate (MTS) reagents. In combi- 
nation with disulfide cross-linking experiments and elec- 
trostatic energy calculations, our data suggest that ions do 
not use the central pathway, but likely permeate through 
the three lateral fenestrations located between the ex- 
tracellular region and the transmembrane pore. 



MATERIALS AND METHODS 

Poisson-Boltzmann calculations 

The zfP2X4 crystal structure (Protein Data Bank accession no. 3H9V) 
was modeled from residue .32 to 361, with missing side chains 
built using CHARMM. The rP2X2 receptor and the human 
P2X5 (hP2X5) receptor homology models were created based 
on the zfP2X4 receptor crystal structure using MODELLER 
software (Eswar et al., 2006). The electrostatic free energy of ion 



interaction along the two pathways was determined by calculating 
AAGi„t = AGion+ehaimei-AG.hannei-AGion. with each Component ob- 
tained from a distinct numerical solution of the finite-difference 
Poisson-Boltzmann equation for the given system. This procedure 
ensures that the contribution from the static field (arising from the 
protein charges) and the reaction field (arising from the spatial vari- 
ations in the dielectric constant) are both incorporated into the 
calculated free energy. The calculations were performed following 
the methods described in Robertson et al. (2008) , using the PBEQ 
module in CHARMM (version .S.5al ) and the set of optimized atomic 
Born radii for amino acids (Nina et al., 1997; Im et al., 1998; 
Brooks et al., 2009). The zfP2X4, rP2X2, and hP2X5 systems 
were aligned with the central pore along the z axis. For the finite- 
difference Poisson-Boltzmann calculations, each system con- 
tained 1 60 X 1 60 X 290 grid points with cell size of 1.0 A focused 
to 0.5 A. The membrane was defined as a 25-A-thick slab centered 
at Z = 0 A with dielectric constant Em = 2, whereas protein and sol- 
vent dielectric constants were set to ep = 2 and es = 80, respectively. 
Ions were treated explicitly inside a box defined from — 15 A < X 
< 15 A, -15A<Y< 15 A, and -20A<Z<90A, and Debye-Hiickel 
ionic screening corresponding to 150 mM of salt was considered 
outside this region. All residues were set to their native proton- 
ation states, with histidines neutral. For the lateral fenestrations, 
AAGint was calculated for accessible points, those positions >2 A 
from a protein atom, within a 20-A radius of Z = 25 A. The artifi- 
cial opening of the central pore was constructed by carrying out 
an adopted basis Newton-Raphson minimization with a cylindri- 
cal potential resistant along the central axis of the protein, succes- 
sively increasing the radius by 0.1 A to a final radius of 5 A. 

Expression constructs and transfection 

To avoid unwanted modification at native Cys, all mvitants were 
introduced into the rP2X2-.3T, where all three nondisulfide bridge- 
forming Cys (C9, C.S48, and C4.S0) are mutated to threonine 
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Figure 1 . Two potential pathways for extracellular ions to access 
the transmembrane pore in P2X receptors. (A) Surface representa- 
tion of the zfP2X4 receptor model in the closed state. Each color 
represents one of the three subunits. (B) A sagittal section of the 
closed-state model reveals two potential ion access pathways (central 
and lateral pathways; blue dashed arrows). Electrostatic potential 
surface calculated using the APBS tools (Baker et al., 2001) shows 
two acidic vestibules along these potential ion access pathways 
(upper and central vestibules). The surface is colored based on the 
potential, contoured from —6 kcal mol"' (red) to +6 kcal mol"' 
(blue). Wliite denotes 0 kcal moP'. To complement the electro- 
static contributions of the missing side chains in the crystal struc- 
ture because of weak electron density, the side chains were added 
back to the model, and their geometry was optimized using the 
MODELLER software (Eswar et al. , 2006) . Structural representations 
in all figures were generated using PyMOL (Sclirodinger, LLC). 
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(Li et al., 2008). The parent rP2X2-3T was subcloned into tlie 
pIRES-EGFP RK6 vector (provided by M. Mayer, National Insti- 
tutes of Health, Bethesda, MD) using BamHI and Hindlll sites, 
modified to either single or double Cys mutants using the Quick- 
Change technique (Agilent Technologies) , and confirmed by DNA 
sequencing. For whole cell patch clamp recordings, HEK293 cells 
were split on glass coverslips in a six-well plate (Corning) at 5-10 x 
10* cells/well, cultured in DMEM (Invitrogen) for 10-15 h at 
37°C, transfected with 300-800 ng DNA per well using Fugene 6 
reagent (Roche), and used 16-32 h after ttansfection. For Western 
blotting, ~90% confluent HEK293 cells in six-well pates were trans- 
fected with 4 Jig of DNA constructs/ well using Lipofectamine 2000 
reagent (Invitrogen) and used 48 h after transfection. 

Whole cell patch clamp recording 

Membrane currents were recorded from HEK293 cells using the 
whole cell patch clamp configuration. Membrane voltage was 
clamped at —60 mV with an Axopatch 200B patch clamp ampli- 
fier (Axon Instruments) , currents were filtered at 2 kHz (eight-pole 
Bessel; model 900BT; Freqviency Devices) and digitized online using 
a Digidata 1 32 1 A interface board and pCLAMP 1 0 software (Axon 
Instruments). The standard extracellular solution contained: 
140 mM NaCl, 5.4 mM KCl, 2 mM CaClj, 0.5 mM MgClj, 10 mM 



HEPES, and 10 mM D-glucose, adjusted to pH 7.3 with NaOH. 
The pipette solution contained: 140 mM NaCl, 10 mM EGTA, and 
10 mM HEPES, adjusted to pH 7.0 with NaOH. Bath and ground 
chambers were connected by an agar bridge containing 2 M KCl. 
Extracellular solutions were rapidly exchanged (~50 ms) to ATP- 
containing solutions with and without MTS reagents using a com- 
puter-controlled perfusion system (RSC-200; Bio-Logic). Stock 
MTS reagents (bromide salt; Toronto Research Chemicals) were 
freshly made in HjO at 100 mM (prepared every 4 h) , stored on ice, 
and diluted to the desired concentrations within 2 min of their 
use in experiments. To expedite the screening process, we deter- 
mined the titratable range of the concentration-response rela- 
tion for ATP for each single and double Cys mutant by applying 
three different concentrations of ATP (e.g., 1, 10, and 100 pM). 
For each mutant, MTS accessibility was assessed using an ATP con- 
centration near the EC50, and all the recordings were performed 
for three or more cells. 

Detection of intersubunit disulfide bridge formation 

To obtain a high signal to noise ratio in Western blots, an octa-histi- 
dine tag (Hisg) was attached to each of the double Cys mutants using 
PCR and standard molecular biology techniques. AU such constructs 
behaved similarly to the ones without Hisg in whole cell patch clamp 
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Figure 2. Electrostatic interaction free energy of ions along the central pathways in zfP2X4 and rP2X2. Results for zfP2X4 are shown in 
A and B, whereas those for rP2X2 are shown in C and D. (A and C) AAGint of Na* (black) , Ca^* (green) , and CV (blue) , calculated along 
the central pathway, with the transmembrane region depicted with the yellow slab and the center of the membrane at Z = 0 A. The pore 
radius profile is shown in the bar graph on top. (B and D) AAGint of Na* along the central pathway of the original closed model (black), 
and an artificially opened pore widened to 5-A radius from the molecular threefold symmetry axis (red) . The surface representations of 
the original closed (white) and the forced open (red) conformations are shown. 
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recordings (not depicted). HEK293 cells expressing these double 
Cys mvitants were washed with PBS, resuspended in 100 pl/well of 
S buffer (25 mM Tris, pH 7.4, 137 mM NaCl, .3 mM KCl, 20 mM 
w-dodecyl-p-maltoside. Complete Protease Inhibitor Cocktail [Roche] , 
and 1 mM PMSF), rotated for 30 min at 4°C, and spun at 40,000 rpm 
in a TLA 100.3 rotor (Beckman Coulter) for 30 min. The soluble 
fraction (10 pi) was incubated with nonreducing SDS-PAGE sample 
buffer and resolved on a 4—12% Bis-Tris SDS-PAGE gel (Invitrogen). 
The double Cys mutants were detected by Western blotting using 
1/500 dilution of anti-His monoclonal antibody (Thermo Fisher Sci- 
entific), 1/1,000 dilution of anti-mouse secondary antibody 
conjugated with alkaline phosphatase (Bio-Rad Laboratories), and 
Colorimetric AP substrate (Bio-Rad Laboratories) . 

RESULTS 

The objective of this study was to identify the permeation 
pathway for ions through the extracelhilar region of 
P2X receptor channels. We chose to use the rP2X2 chan- 
nel in this study because: (a) it desensitizes relatively 



slowly, thereby facilitating electrophysiological accessi- 
bility studies; (b) ion accessibility of the pore has been 
well characterized (Li et al., 2008, 2010; Kracun et al., 
2010); and (c) it has high sequence identity to zfP2X4 
in the transmembrane and extracellular regions (^^50% 
identity) . To facilitate subsequent computational analy- 
sis and experiments, we also created an rP2X2 homol- 
ogy model based on the closed-state zfP2X4 structure. 

Electrostatic free energy of cations along the potential 
pathways in the closed state 

To assess the energetics of ions along the two pathways, 
we calculated the electrostatic free energy of transferring 
Na^ or Ca'* from aqueous solution to the central pore of 
zfP2X4. The electrostatic free energy profile, shown in 
Fig. 2 A, includes both the static field arising from the 
protein charges, as well as the reaction field resulting 
from the different dielectric regions. In this closed-state 
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Figure 3. Electrostatic interaction free energy of ions along the lateral fenestrations in zfP2X4 and rP2X2. Results for zfP2X4 are shown 
in A, whereas those for rP2X2 are shown in B. The ion-accessible positions are shown as orange spheres near the extracellular membrane 
interface around Z = 15-25 A (left). Free energy plots (AAGim) for cross sections through the lateral fenestrations (Z = 16, 18, and 20 A) 
are shown on the right for Na* (top) and Cr (bottom) . Black squares indicate AAGint < — 2 Kcal mor. 
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structure, the interaction energy of cations is prohibitive 
near the extracellular entrance of the central pathway 
(55 A < Z < 85 A) . In contrast, the center of the extracel- 
lular domain (30 A < Z < 55 A) is extremely favorable for 
both Na* and Ca^^, whereas the energy of CP is repulsive 
at each point along the central axis. Both of these obser- 
vations result from a strongly electronegative vestibule, 
lined by six negative residues (E98 and D99), which cor- 
responds to the Gd^^-binding site in the zfP2X4 struc- 
ture. A simple expanded pore model was created to probe 
the degree of widening needed to reduce the electro- 
static reaction field barrier and assess the static charge 
contribution to ion-interaction energies at the central axis. 
This was done by imposing a repulsive cylindrical poten- 
tial restraint at the center of the pore and carrying out 
energy minimizations of the structure while successively 
increasing the cylinder radius. The resultant model has 
a pore that is 10 A in diameter and maintains the overall 



distribution of residues lining the pore in the closed 
structure. In this model, the electrostatic free energy 
for Na^ is favorable along most of the central pathway 
(Fig. 2 B), making it a possible pathway for cation per- 
meation in zfP2X4. However, when we performed the 
same calculations for the rP2X2 homology model, there 
were still significant electrostatic barriers along the pore 
that occur because of the lack of acidic residues in the 
sequence of this homologue (Fig. 2, C and D) . This sug- 
gests that although the central pathway is a favorable 
option in zfP2X4, the sequence variation in rP2X2 cre- 
ates an unfavorable electrostatic environment that is un- 
likely to support cations along the central axis. 

The electrostatic accessibility of the lateral fenestra- 
tions was also examined in the zfP2X4 structure and 
rP2X2 model (Fig. 3, A and B) . We found that for both 
of these closed structures, there was a clear pathway 
connecting the outside solution to the central pore that 
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Figure 4. Accessibility of MTS reagents along the lateral fenestration. (A) The mutated residues are shown in stick representation 
(blue, carbon atoms; red, oxygen; dark blue, nitrogen) , and the remaining region is shown in cartoon representations (gray) . A trans- 
membrane pore4ining residue, T3.S6, is shown in orange. Representative residues are labeled. (B-E) Accessibility of two different-size 
MTS reagents in the presence of ATP. Superimposed scaled current traces recorded in response to ATP without (black) and with (blue) 
MTSET or MTS-TPAE applications (blue bars) show no modification for Q52C (B) and S326C (C) , and modification for I328C (D) and 
T336C (E) . (F) TR-MTSEA is accessible to I328C in the open state but not in the closed state. Exposing the channels to TR-MTSEA for 
10 s in the closed state did not modify the ATP-induced current, whereas a 6-s exposure in the presence of ATP modified the current. 
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was electrostatically accessible by Na* (<1 kcal moP^). 
Notably, CP was not energetically stable along these path- 
ways. Small openings of the windows would be sufficient 
to remove any small electrostatic barriers created by the 
protein reaction field and provide an energetically 
favorable pathway for cations. These calculations suggest 
that it is electrostatically possible for cations to move via 
the lateral fenestrations, even in the closed-state struc- 
ture of both zfP2X4 and rP2X2. 

Accessibility and structural changes around the 
lateral fenestrations 

The lateral fenestrations consist of two strands connect- 
ing the two membrane-spanning a helices within each 
subunit to the large extracellular ligand-binding region. 
In contrjist to the long and narrow central pathway, these 
lateral fenestrations are relatively wide and shallow in 
the closed-state structure (Fig. 1 B) . We set out to explore 
whether these lateral fenestrations serve as ion permeation 
pathways by substituting Cys residues at each of the 1 7 resi- 
dues that make up the walls of the fenestrations (Fig. 4 A) , 
and tested whether their reaction with MTS reagents has 
detectable effects on ATP-activated currents. In the closed- 
state structure, these residues exhibit extensive solvent 
exposure and therefore should react with MTS reagents. 
However, if the wide and shallow features of the lateral 
fenestrations are maintained in the open state, one might 
imagine that MTS reagents would not have a dramatic 
effect on ion conduction, but might alter gating of the 
channel because these strands connect the ligand-binding 
region to the pore of the channel. For these experiments, 
we externally applied 2-(trimethylammonium)ethyl MTS 
(MTSET; mol wt: 278 D) and 2-(tripentylammonium) ethyl 
MTS (MTS-TPAE; mol wt: 447 D), two different-size re- 
agents that are water soluble, membrane impermeable, 
and positively charged. Because both reagents react rel- 
atively rapidly with T336C within the gate region of 
rP2X2 in the presence of ATP ('^10* M^^s ), and cause 
reductions of ATP-induced currents (Fig. 4 E) (Li et al., 
2008, 2010), the extracellular pathway to the pore must 
be wide enough for these MTS reagents to access the 
ion permeation pathway. 

We examined the effects of these two MTS reagents on 
ATP-activated currents for each of the 15 Cys mutants that 
retained functional activity when expressed in HEK cells. 
Y55C and Q,56C failed to form functional channels and 
were not studied further Whole cell patch clamp record- 
ings from these mutants revealed that the application of 
either MTSET or MTS-TPAE has negligible effects for 
most Cys mutants (Figs. 4, B and C, and 5 B) . Modest cur- 
rent reduction was observed with MTSET for the I328C 
mutant (Fig. 4 D, left traces), whereas small current po- 
tentiation was recorded with MTS-TPAE (Fig. 4 D, right 
traces) , suggesting that 1328 in the lateral fenestration 
is accessible to both MTS reagents. The rate of MTSET 
modification at I328C was --5 x 10^ M^^s"^ (n = 8), 



slower than the rates observed for MTSET modification 
of most pore-lining residues within TM2. Consistent with 
our expectations from the zfP2X4 structure, these results 
demonstrate that only one residue forming the lateral fen- 
estrations exhibits detectable reaction with MTS reagents. 

Although the zfP2X4 structure shows that the resi- 
dues we tested in the lateral fenestrations are solvent 
exposed, and therefore should react with MTS reagents, 
an alternate possibility is that these lateral fenestrations 
collapse when the channel opens and their solvent acces- 
sibility decreases. To test this possibility, we examined 
whether the larger mushroom-shaped Texas red-2-sulfon- 
amidoethyl MTS (TR-MTSEA; molwt: 744D) can react with 
I328C in either the absence or presence of ATP. Although 
TR-MTSEA does not react with I328C in the absence 
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Figure 5. Quantitative comparisons of current modifications by 
MTS reagents along the lateral window. (A) Superimposed cur- 
rent traces of T336C recorded in response to ATP without (black) 
and with (blue) MTSET application (blue bar) using a 2-min in- 
terval between traces to illustrate how modification is quantified. 
Traces were scaled using the current amplitude 1.5 s after ATP 
application. Percent modification was calculated according to: 
modification (%) = [(b/a) — 1] x 100, where b is the normalized 
current after a 6-s MTSET application, and a is the current at the 
corresponding time in the control trace. (B) Averaged current 
modifications by MTSET (left) or MTS-TPAE (right) at each Cys 
mutant (ra = 3-9) . Error bars are SEM. Orange bars indicate posi- 
tions where current reduction or potentiation was observed. 
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of ATP, consistent with its large size (^6 x 13 x 16 A) 
and the dimensions of the lateral fenestrations in the 
closed-state zfP2X4 structure, it reacts rapidly in the pres- 
ence of ATP (Fig. 4 F) . Interestingly, the reagent initiaJly 
potentiates ATP-activated current, similar to MTS-TPAE, 
but then causes complete inhibition of ATP-activated cur- 
rents. One possible explanation is that the reaction with 
each subunit stabilizes the open state of the pore while oc- 
cluding a lateral fenestration, and that complete inhibi- 
tion is not detected until all three fenestrations are blocked. 
Regardless of the precise mechanisms, these results sug- 
gest that the lateral fenestrations become larger upon ATP 
binding; therefore, it is a likely pathway for ions to enter 
and exit the pore. 

Accessibility and structural changes in the central pathway 

To explore whether the central pathway might serve as 
a permeation pathway, we introduced Cys residues at 26 
positions (Fig. 6 A) corresponding to residues at the cen- 
tral threefold axis in the zfP2X4 structure and tested for 
modification by MTS reagents. For 24 of these positions, 
both MTSET and MTS-TPAE had no effect on ATP-acti- 
vated currents (Figs. 6, B-E, and 7 C) . In the acidic central 
vestibule along the central pathway, two mutants (I317C 
and H319C) were modified by MTSET, but not by the 
larger MTS-TPAE. For I317C, MTSET increased the rate 
of deactivation (Fig. 7 A, middle traces), whereas for 
H319C, MTSET potentiated the ATP-induced current 
by ^40% and slowed the rate of current deactivation 
(Fig. 7 B, middle traces) . In the caise of I317C, subsequent 



application of ATP 3 min after the removal of MTSET 
resulted in ATP-activated currents that retained rapid 
deactivation, suggesting that the initial effects of MTSET 
are the result of covalent modification of 13 1 7C (Fig. 7 A, 
right traces). In the caise of H319C, the potentiation of 
ATP-activated currents by MTSET did not recover cifter 
the removal of the MTS reagent, and a subsequent appli- 
cation of MTSET in the presence of ATP wais without ef- 
fect (Fig. 7B, righttraces) , indicating covalentmodification 
of H319C by MTSET. The rate of MTSET modification at 
H319C was ^3.5 x 10^ M^^s"^ (n = 5), consistent with this 
residue being solvent exposed in the presence of ATP. 
In both cases, preapplication of MTS-TPAE did not af- 
fect the MTSET-dependent current modification (Fig. 7, A 
and B, left and middle traces) , suggesting that the larger 
MTS-TPAE cannot access these residues. 

To explore whether the central pathway undergoes a 
large conformational change upon ATP binding to cre- 
ate a wide route for extracellular ions to permeate, we 
engineered five pairs of double Cys mutants along the 
central pathway (Fig. 8 A) , which based on the closed-state 
zfP2X4 structure would be predicted to form intersub- 
unit disulfide bonds. Because P2X receptors are trimers, 
each Cys pair could give rise to a total of three disul- 
fides. Such bridges should prevent a large motion because 
of conformational restraints caused by introducing in- 
tersubunit disulfide bridges. These double Cys mutants 
were expressed in HEK cells, and their ability to form 
disulfide bridges was assessed by SDS-PAGE in the ab- 
sence of a reducing agent, followed by Western blotting. 
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Under mild air-oxidizing conditions, the band for the 
wild-type rP2X2 migrated at the appropriate molecular 
weight for the monomeric species (^^53 kD), indicating 
that there is no intersubunit disulfide bridge formation 
in the parent construct (Fig. 8 B, lane WT) . On the other 
hand, all double Cys mutants formed SDS-resistant trimers 
(^^150 kD) or dimers (^100 kD), suggesting that inter- 
subunit disulfides are formed, at least to some extent, in 
all the mutants (Fig. 8 B, lanes I-V) . 

Under identical air-oxidizing conditions, we performed 
whole cell patch clamp recordings from HEK cells express- 
ing each pair of double Cys mutants. We chose to apply 
ATP (10 pM for the pairs I and III-V, and 30 pM for the 
pair II) for relatively short durations (Is) at 2-min inter- 
vals because under these conditions, we observed rela- 
tively litde rundown upon repeated applications of ATP 
(Fig. 9, left and middle traces) . Notably, when cells were 
perfused with 10 mM DTT for 5 min to reduce disul- 
fides, the top pair along the central pathway (Fig. 9, I) 
showed substantial current inhibition (^70%) , whereas 
currents in three other pairs (II-IV) were unchanged. 
The data for bridge I indicate that the channel activity 
detected for the outer pair of Cys mutants under oxidiz- 
ing conditions must be mediated, at least in part, by the 
Cys mutants whose subunits were covalently cross-linked 
by intersubunit disulfides. Thus, we conclude that restrict- 
ing the subunit-subunit movements at the top of the cen- 
tral pathway is compatible with ion permeation through 
the rP2X2 channel, supporting the idea that a large open- 
ing of the central pathway is unlikely. The interpreta- 
tion of the three bridges where DTT has no effect (Fig. 9, 
II-IV) is ambiguous because, as pointed out in the previ- 
ous studies on channels (Kobertz et al., 2000), it is un- 
clear whether the currents we observed are mediated by 
the population of disulfide cross-linked channels. We tried 
to drive disulfide bridge formation to completion by treat- 
ing HEK cells with 0.3% H2O2 for 10 min, but no signifi- 
cant increase in disulfide formation was obtained (not 
depicted) . In the case of the bottom pair (V) , we observed 
large current potentiation when the cell was exposed to 
DTT, implying that the disulfide bond impairs a confor- 
mational change involved in channel opening. This 
bridge is actually positioned at the top of the lateral fenes- 
tration, and these results are therefore consistent with 
experiments with TR-MTSEA, which suggest that the lat- 
eral fenestrations become larger when the channel opens. 

DISCUSSION 

In this study, we explored the pathway for ions to per- 
meate through the extracellular region of P2X receptor 
channels. The zfP2X4 crystal structure in the closed state 
reveals the presence of two potential pathways. Our 
Poisson-Boltzmann calculations suggest that substantial 
structural changes (>5 A from the central axis) would 
be required for ions to permeate through the central 
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pathway (Fig. 2) . Although our accessibility experiments 
show that ions up to the size of MTSET can reach the 
central vestibule (Fig. 7) , we failed to observe any modi- 
fication of residues in the rest of the central pathway 
(Fig. 6) . In addition, large MTS reagents like MTS-TPAE 
can readily access the transmembrane pore but cannot ac- 
cess either the central or upper vestibules (Figs. 6 and 7) , 
suggesting that the central pathway is unlikely to be an 
ion permeation pathway. Importantly, engineered di- 
sulfide bridges that should constrain movement at 
the upper end of the central pathway do not significantly 
perturb channel function (Figs. 8 and 9), excluding the 
possibility that the central pathway undergoes a large 
dilation upon activation. In contrast, our computational 
analysis shows no significant electrostatic barriers for 
cations to enter from the lateral fenestrations (Fig. 3) , 
and our accessibility results identify one position where 
both MTSET and MTS-TPAE can access an introduced 
Cys in the presence of ATP (Fig. 4). 

Enlargement of the lateral fenestrations in the open state 

One expectation of accessibility studies with ion chan- 
nels is that the introduction of a positive charge with an 
MTS reagent should reduce macroscopic current if the 
Cys is positioned along a relatively narrow permeation 
pathway. These expectations were met in previous stud- 
ies identifying pore-lining residues and the gate within 
TM2 of the rP2X2 receptor channel (Li et al., 2008, 2010) . 
The relative scarcity of reactive positions within the cen- 
tral pathway in the present study is not surprising if it 
does not contribute to the ion permeation pathway. The 
central pathway is '^70-A long and includes extended 
regions where packing within the protein minimizes its 
diameter, and at least at its outer end does not undergo 
a substantial structural change during gating, as inferred 
from the engineered disulfide that retains activity. How- 
ever, after screening 15 Cys mutants within the lateral 
fenestrations with two different-size MTS reagents, we 
found only one mutant (I328C) that decreased ATP- 
induced current upon MTSET application (Fig. 4) , and 
in this instance, the rate of modification was relatively 
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Figure 9. Disulfide bridge formation in the central pathway pre- 
serves channel activity. Two short pulses of ATP applications (1 s; 
black arrows) with 2-min intervals showed inward currents with- 
out significant rundown. External DTT application (5-min dura- 
tion) reduced the current amplitude of the mutant pair I and 
potentiated mutant pair V (red asterisks) . 
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moderate (^5 x 10^ M^^s"^). Although the modest 
modification rate could be the consequence of various 
local factors (e.g., dynamical, structural, or energetic) 
suppressing the reactive thiolate ion in this region, the 
lack of effect at the other 14 positions within the lateral 
fenestration suggests that the diameter of the opening 
is actually quite large in the open state. In this case, 
many of the fiinctionally silent positions may have reacted 
with MTS reagents, but the size of the opening prevents 
the positive charge from attenuating ion throughput. 
This interpretation would be consistent with a recent 
MTS accessibility study on pannexin hemichaiinels (Wang 
and Dahl, 2010), which failed to identify any reactive 
positions within the four TM helices, even though the 
recent crystal structure of a potential structural analogue, 
the gap junction channel (Maeda et al., 2009), would 
predict many to be highly solvent exposed. Enlargement 
of the lateral fenestrations when P2X receptor channels 
open is strongly suggested by accessibility experiments 
with large MTS reagents like MTS-TPAE (--6 x 12 x 13 A) 
and Texas red (^6 x 13 x 16 A), which can rapidly enter 
the pore and react with pore-lining residues (Fig. 4) (Li 
et al., 2010), requiring that the fenestrations be signifi- 
cantly larger than these reagents. Indeed, our present 
experiment with TR-MTSEA and position I328C within 
the lateral fenestration suggests that it becomes larger 
as the channel opens. The breathing of the lateral fen- 
estrations with opening would also be consistent with 
the effects of the intersubunit disulfide bridge in this 
region (V) . Reduction of this bridge with DTT enhances 
macroscopic current (Figs. 8 and 9), suggesting that 
constraining motions of the lateral fenestration prevents 
the opening of the channel. 



Ion selection In the lateral fenestrations 
In identifying the lateral fenestration as the permeation 
pathway for ions to enter or exit the pore, it becomes 
interesting to consider whether this pathway may con- 
tribute to the ion selectivity of P2X receptor channels. 
Our Poisson-Boltzmaiin calculations reveal that the elec- 
trostatic interaction energy along this short (^15-A) path- 
way is relatively favorable for Na^ permeation ( AAGim = 
-1 to+lkcal mol *), butnotfor Cr (AAGim > +10 kcal 
moP^) (Fig. 3), suggesting that the overall negative po- 
tential in this area may contribute to selecting cations 
over anions. Interestingly, the P2X5 receptor has been 
reported to also be permeable to anions (pCl /pNa^ = 
--0.5; Ruppelt et al., 2001; Bo et al., 2003), and the elec- 
trostatic potential in a hP2X5 receptor homology model 
is significandy positive around the lateral fenestration 
(^^+10 kcal mol Fig. 10), supporting the hypothesis 
that the electrostatics along the lateral fenestration plays a 
role in ion selection. 

Role of the acidic central vestibule 

Why do P2X receptors contain a large acidic vestibule 
in the middle of the central pathway? In our MTS ac- 
cessibility studies, we found that the smaller MTSET 
(^6 X 6 X 10 A) can access the acidic central vestibule 
through the lateral fenestrations and potentiates currents 
(H319C) or accelerates the rate of deactivation (I317C) 
(Fig. 7, A and B). These data suggest that extracellular 
ions can reach the central vestibule and allosterically 
modulate conformational changes in the extracellular 
domain involved in opening the pore. Supporting this, 
a trivalent cation (Gd^^) has been shown to bind in the 
central vestibule of zfP2X4 and inhibit ATP-induced 
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Figure 1 0. Electrostatic potential in the lateral fenestrations may contribute to selecting cations over anions. Sagittal sections of the 
electrostatic potential surface calculated using the APBS tools (Baker et al., 2001). (A) zfP2X4 with supplemented side chains. Homol- 
ogy models of rP2X2 (B) and hP2X5 (C). Surface is colored based on the potential, contoured from —6 kcal mol"' (red) to +6 kcal 
mol~' (blue). White denotes 0 kcal mol"'. 
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currents (Kawate et al., 2009) . It is also tempting to spec- 
ulate that the acidic surface may serve a role in concentrat- 
ing cations, facilitating their influx to the transmembrane 
pore. Interestingly, previous studies have shown that in 
large-conductance Ca'^-activated (BK) and KcsA 
channels, intracellular is locally concentrated by a 
ring of negatively charged glutamate residues at the 
entrance of the intracellular vestibule (Brelidze et al., 
2003; Nimigean et al., 2003). In the study on BK chan- 
nels, the local concentration of at the internal mouth 
of the pore is estimated (based on the reversal poten- 
tials) to increase from 150 to 500 mM as a result of the 
negative charges. If the central vestibule in the extracel- 
lular region of a P2X receptor increases the local con- 
centration of Na^ in a similar manner to BK channels, 
this acidic central vestibule should have a profound 
effect on cation/ anion selectivity. 
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